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Abbreviations 

Some abbreviations are used in this chapter (by alphabetic order): 
bct=body-centered-tetragonal, 

BS=Breakdown of Symmetry, d=dimension, fcc=face-centered cubic lat- 
tice, FSS=Finite Size Scaling, GS=Ground State, hcp=hexagonal-close- 
packed lattice, KT=Kostcrlitz-Thouless, L=system size, MC=Monte Carlo, 
MCRG=Monte Carlo Renormalization Group, N— number of components 
of the spin S, NN=Nearest Neighbors, NNN=Next Nearest Neighbors, 
RG=Renormalization Group, STA=Stacked Triangular Antiferromagnctic 
lattices, STAR=Stacked Triangular Antiferromagnetic lattices with Rigid- 
ity, T=Temperature, T c =critical Temperature, V/v,p=Sticfcl model. 

1. Introduction 

We present in this chapter a review on recent numerical studies dealing 
with frustrated vector spin systems in two and three dimensions. A system 
of spins is frustrated when all interactions between spin pairs cannot have 
simultaneously their optimal values. In other words a system is frustrated 
when the global order is incompatible with the local one, a definition ap- 
plicable in a broader sense and not restricted to spins. For spin systems 
frustration has the consequence that the ordered state is different from the 
collinear order found for common unfrustrated antiferromagnets or ferro- 
magnets. To study and to classify the phase transitions between the ordered 
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2 D. Loison 

and less ordered states the symmetries these states play an important role. 
The critical behavior of second order transitions is principally governed by 
the change of symmetries. Therefore even very different systems like super- 
fluid helium and XY spin systems can have the same critical properties. 
This fact is called universality and our objective is to analyze the corre- 
sponding universality classes for frustrated spin systems. We will mainly 
review magnetic models since they are the easiest to analyze numerically 
and theoretically. Nevertheless the results should be valid for any system 
belonging to the same universality class. 

During the last decade important progress has been made in the under- 
standing of the physics of frustrated spin systems. For example there is now 
convincing evidence on the genuine first order nature of the phase transition 
for XY and Heisenberg spins for the three dimensional Stacked Triangu- 
lar Antiferromagnet (STA). This still contradicts the latest renormalization 
group expansion based on resummation. We think that studying phase tran- 
sitions of vector spins theoretically, numerically and experimentally should 
have implications beyond this special field, that is for the understanding 
and the theory of phase transition in general. 

We will concentrate our attention to the phase transition of the physical 
XY and Heisenberg spin systems in two and three dimensions. However, 
to understand these systems we have to analyze also the phase transitions 
of frustrated spins of N components, where N takes all integer values from 
1 to oo and not only 2 and 3 for XY and Heisenberg spins. We will also 
present studies for dimension d varying between two and four. Then we 
will review the particular case of strictly two dimensions where topological 
defects have a dominant role. 

Frustrated Ising models are reviewed also in this book by Diep and 
Giacomini (chapter 1) and by Nagai, Horiguchi and Miyashita (chapter 2). 

Since most of numerical simulations presented here use the Monte Carlo 
method, a short appendix @ at the end is devoted to this technique. In ad- 
dition, since the renormalization group is fundamental to the understanding 
of phase transitions, a small appendix © is added to discuss the methods 
used here for frustrated spin systems. 

2. Breakdown of Symmetry 

We first briefly review the fundamental concept of the reduction or the 
Breakdown of Symmetry (BS) in the transition from the high- to the low- 
temperature phases. The classification of phase transitions in universality 
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classes is based on this concept. 



2.1. Symmetry in the high-temperature region 



We will consider the Hamiltonian: 



m 



(i) 



where Si are N component classical vectors of unit length and the sum 
is usually restricted to the nearest neighbors or at least to short range 
interactions. The symmetry of this Hamiltonian is O(N) or equivalently 
Zi (g) SO(N) where O(N) is the orthogonal transformation group of TV- 
dimensional Euclidean space. For SO(N) the determinant is unity, and the 
Ising symmetry Z2 corresponds to the mirror operation. If we add other 
terms, like dipolar interactions or anisotropies, the symmetry group will be 
reduced. Some other changes, like a cubic term (S^.S^) 3 , will not reduce 
the symmetry. In many experiments anisotropies will reduce Hciscnbcrg 
(N = 3) symmetry to XY (N = 2) or Ising (N — 1) symmetry. Long-range 
interactions could also be present and the interpretation of experimental 
results could be problematic (see the crossover section inEJ. 

We will also encounter the Potts symmetry Z q . The Potts model has 
the Hamiltonian: 



5 qiqj refers to the q states Potts spin with S qiQj = when qi ^ qj and 

SqiU = 1 wnen ft = Qj- 

In addition we have to consider the symmetry of the lattice which is 

the sum of all symmetry elements which let the lattice invariant. For ex- 
ample, the triangular lattice has C% v symmetry and the square lattice C4 
symmetryff] Therefore we have to take into account the total BS composed 
of the 0(N) spin rotations and the symmetries of the lattice. 

2.2. Breakdown of symmetry for ferromagnetic systems 

If Ji in is taken positive the Ground State (GS), for any lattice, is 
ferromagnetic with collinear spins pointing in the same direction. If J\ is 
negative (antiferromagnetic) but there is no frustration like for the square 
lattice, the GS is collinear with alternating direction. Then the BS will be 
identical to the ferromagnetic case. The symmetry at low temperature is 




(2) 



(ij) 
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0(N — 1) for the spins and the lattice symmetries are preserved. Therefore 
the BS is 0(N) — > 0(N — 1) where O(N) is the symmetry of the phase 
at high temperatures and 0(N — 1) the symmetry at low temperatures. 
We will denote this BS 0{N)/0{N - 1). Keeping in mind that O(0) = 1, 
0(1) = Z 2 , SO(l) = 1 and O(N) = Z 2 ® SO(N) for N > 2, the BS 
0(N)/0(N - 1) can be written as Z 2 for JV = 1, 50(2) for N = 2, and 
SO(A)/S*0(A - 1) for N > 2. 

All the ferromagnetic systems (cubic, stacked triangular, . . . ) with a 
Hamiltonian of type Q have an identical BS and consequently belong to 
the same universality class. The only relevant variables are the dimension 
of the spin space, i.e. the number of components of the spin N, and the 
dimension of the real space d. 

For Potts ferromagnetic system with a Hamiltonian of the type © and 
J\ > 0, the Potts symmetry is broken in the low temperatures phase and 
the BS is simply Z q . Again this does not depend on details of the system. 
This model has a second order phase transition for q < q c and a first order 
phase transition for q > q c . We know that q c = 4 in two dimensions and 
q c < 3 in three dimensions. Therefore, in three dimensions, for any q > 3 
the transition will be of first order. 

In Table ^ we g ive the value of the critical exponents calculated by the 
Renormalization Group^ (RG) depending on N and q in three dimensions 
(d = 3). Even if > 3 does not correspond to physical systems, we will 
get very useful information for frustrated systems. 

Table 1. Critical exponents for the ferromagnetic systems in three di- 
mensions calculated by RG. ( a 'We cannot define exponents in a first-order 
transition, however in the case of a weak first-order transition the exponents 
found with MC and in experiments must tend to these values, '''calculated 

by y/v = 2~ri. 



BS 


a 


P 


7 


V 


V 


O(l)/O(0) = z 2 


0.107 


0.327 


1.239 


0.631 


0.038 


0(2)/0(l) = SO(2) 


-0.010 


0.348 


1.315 


0.670 


0.039 


0(3)/0(2) = SO(3)/SO(2) 


-0.117 


0.366 


1.386 


0.706 


0.038 


0(4)/0(3) = SO(4)/SO(3) 


-0.213 


0.382 


1.449 


0.738 


0.036 


0(6)/0(5) ee SO(6)/SO(5) 


-0.370 


0.407 


1.556 


0.790 


0.031 


1st order< a > (Z q if q > 3) 


1 





1 


1/3 


-l(«0 
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2.3. Breakdown of symmetry for frustrated systems 

The BS in frustrated systems is more complicated. We will present examples 
where the 0(N) symmetry is reduced to 0(N - 1), 0(N - 2), • ■ ■ , 0(0) 
(completely broken in the last case). In addition it is possible that the 
lattice symmetry is also reduced, usually giving an additional Z q (mostly 
<7 = 2 and q = 3) broken symmetry. Therefore the possible BS are Siattice ® 
0{N)/0{N - P) with Siattice = 1, Z 2 or Z 3 and P varies from 1 to N. 

2.3.1. Stacked triangular antiferromagnetic lattices 

We consider here the stacked triangular antiferromagnetic lattice (STA) 
with nearest-neighbor (NN) and next-nearest-neighbor (NNN) J 2 antiferro- 
magnetic interactions, J\ and J%, respectively. There is no frustration along 
the z axis (the direction of stacking) and therefore the interaction can be 
ferromagnetic or antiferromagnetic along z. We will explain this case in 
detail since it shows many phenomena appearing in frustrated systems. 

• J2 = 0: Consider the case without second neighbor interaction 
(J2 = 0). It is not possible for all three spins at the corners of 
a triangle to have the optimal antiparallel orientation which would 
minimize the energy of individual pair interactions. The resulting 
compromise in the case of vector spins is the so-called 120° struc- 
ture, as shown in Fig. Q](see chapter 1). In frustrated systems, local 
minimization of the energy is not compatible with the global energy 



minimum (or minima) . A formal definition proposed by Toulouse"' 
(and also VillairPI) in the study of spin glasses states that a geome- 
try is frustrated if the sign of the product of exchange interactions 
Ji around a plaquette C 



is negative (where Ji < implies antiferromagnetic interactions). 
An antiferromagnetic triangular plaquette is thus frustrated as 
it involves a product of three Ji < 0. The triangular lattice is 
fully frustrated since all plaquettes satisfy this rule. The princi- 
pal effect of the frustration here is that it gives rise to a non- 
collinear magnetic order. This spin-order GS (I) is planar and sta- 
ble as long as next nearest neighbor interaction is small, that is 
< J2/J1 < 0.125. 




$ c = sign[]J Ji] 



(3) 



ieC 
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Fig. 1. Ground State I of the antiferromagnetic triangular lattice (first neighbors only). 
The GS has a 120° spin structure. 



The BS between the disordered phase at high temperatures and the 
phase at low temperatures is 0(N)/0(N — 2) which is equivalent 
to Z 2 ® 50(2) for N = 2 and SO(N)/SO(N - 2) for N > 2. We 
notice that the exact value of the angle between the spin directions 
plays no role for the phase transition. 

This syst em has been extensively stud ied in three dimensions for 
XY spin a I I I I I for Heisenberg spin p l " ! u l ' ' and for the un- 
physical systems N = 6 in Ref. [13] and N = 8 in Ref. [8]. The two 
dimensional systems have also b een the subject of several studies 
for XY spin J 1U7 l 1U8 l 1U9 l 11U l 111 land also for Heisenberg spinsPS 
0.125 < J2/J1 < 1: Consider the presence of an antiferromagnetic 
second neighbor interactions (J2) in the xy plane. The GS can be 
determined by minimizing the energy after a Fourier transform!^ 
For < J2/J1 < 0.125, the GS is still the 120° structure. For 
0.125 < J2/J1 < 1 the GS is degenerate with 9 taking any value 
between and 2ir (see Fig. |2J). However the degeneracy will be 
lifted by thermal fluct uations (spin waves) and only a collinear GS 
(II) will be chosen.^lThis phenomenon, called "order by disorder" 
following Villainp3 is general in frustrated systems and we will see 
other examples (fee, hep) later. There are three ways to choose the 
parallel spins and the BS of the lattice symmetry C% v is a three- 
state Potts symmetry Z 3 . The total BS is Z 3 ® 0(N)/0(N - 1). 
For XY spins (N = 2) it is equivalent to Z 3 ® 50(2), and to 
Z 3 ® SO(N)/SO(N - 1) for N > 2. 

Numerical studies have been d one fo r the three-dimensional case 
by Loison, Diep and BoubcheurP^^ 

Jil J\ > 1: Now consider that J2/J1 > 1. The ground state is also 
degenerate, but this degeneracy is lifted by thermal fluctuations. 
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Fig. 2. Ground State (II) of the antiferromagnetic triangular lattice with first (Ji) and 
second ( J2)-neighbor antiferromagnetic interactions. 0.125 < J2/J1 < 1. (a) 8 can take 
any value between and 2ir. (b) at T > only the "most collinear" GS is chosen. 



The GS (III) is no more collinear, but still planar (see Fig. |3J) and 
there is still three ways to choose the parallel spins, a is determined 
bycos(a) = — 0.5-(l + J1/J2) and can be incommensurable, i.e. does 
not correspond to a rational value. Furthermore a varies slowly as a 
function of the temperature. The BS is Z 3 ®0{N) fO(N - 2) which 
is equivalent to Z 3 ® Z 2 <8 50(2) for N = 2 (XY spins), to Z 3 ® 
SO(3) for TV = 3 (Heisenberg spins), and to Z 3 ®SO(N)/SO(N-2) 
if N > 4. Numerically the incommensurable angle is problematic. 
Indeed we have to impose periodic boundary conditions and since 
a varies with the temperature, the size chosen compatible with 
a(T = 0) is no more compatible with a(T > 0) < a(T — 0) 
for higher temperatures. Therefore the use of Finite Size Scaling 
technique (FSS, seelU) to calculate the critical exponents will be 
problematic!^ In addition the GS III may not be stable under the 
new boundary constraint. Then the Potts symmetry could not be 
broken and the BS could be just 0(N)/0(N — 2), and does not 
belong to the same universality class as previously. We will have an 
akin problem for helimagnets (see later) and for triangular lattices 
with two distinct nearest-neighbor interactions! ^^ 
This model has been studied in three dimensions by Loison, Diep 
and Boubcheur! 7 I H | 

Other BS: We have seen the three BS which appear between the 
disordered phase and the GS. In addition two other transitions 
appear for a small range of J2/J1 near 0.125 and 1: between the 
GS I and II for 0.120 < J 2 / Ji< 0.125, and between the GS II 
and III for 1 < J2/J1 < 1.05 M ' ' This is a general scheme: due 
to thermal fluctuation the most collinear state is favored when the 
temperature increases. Looking at the symmetries we can get the 
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Fig. 3. (a) Ground State III of the antiferromagnetic triangular lattice with first (Ji) 
and second (J2) -neighbor antiferromagnetic interactions. 1 < J2/ ' J\. (b) Phase diagram 
in the (J2/ Ji,T) space. 

order and even the universality class of these transitions. 
The GS III and II are compatible in the sense that at the transition 
a will go smoothly to zero and the transition could be of second 
order. The BS between III and II is 0(N - 1)/0(N - 2) and con- 
sequently the transition should be N — 1 ferromagnetic type. 
The GS I and II are incompatible because there is no way to go 
smoothly from one to the other and the transition should be of first 
order. 

2.3.2. bet Helimagnets 

We consider the body-centered-tetragonal (bet) helimagnets. The compe- 
tition between the J\ and J2 interactions gives rise to a helical ordering 
along the z axis (see Fig. This GS is characterized by the turn angle 
a between spins belonging to two adjacent planes perpendicular to the z 
axis, a is given by the formula cos(a) = — J2/J1 and a decreases slowly as 
function of the temperature. There is no breakdown of the symmetry of the 
lattice, the GS is non collinear but planar and the BS is 0(N)/0(N — 2), 
i.e. identical as the STA with first nearest neighbor interaction only. 

Numerically we have a similar problem as for the STA with large NNN 
interaction: a(T) varies as a function of the temperature, but with periodic 
boundary conditions a constraint is present. 

Nevertheless effect of this constraint should not be too strong. Indeed a 
small constraint will not break a symmetry of the lattice and for the XY, 
it will not change the BS of the spin rotation because all the symmetries of 
the rotation group are already broken. For Heisenberg spins the GS could 
become non coplanar, but this is unlikely. Therefore we should find the 
same universality as that of the STA. The only difference will be a new 
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correction to the scaling laws for a second order phase transition. But for 
a first order one the boundary conditions would no t matt er. 



Fig. 4. a. Interactions for the bet. b. Ground State for the bet for J2/J1 = 0.5. 



2.3.3. Stacked Ji~J 2 square lattices 

The J\-Ji simple cubic lattice is made by stacking along the z axis the 
square lattices with antiferromagnetic (or ferromagnetic) first nearest neigh- 
bor interaction (Ji) and antiferromagnetic second nearest neighbor inter- 
action (J2) in the xy plane. For J2/J1 > 0.5 the GS is degenerate but 
by "order by disorder" only collinear configurations appear (see Fig. [SJ- 
There are two ways to place the parallel spins (following the x or y axis) 
and a Z 2 Ising symmetry of the lattice symmetry is broken. The BS will 
be Z 2 * 0(N)/0(N - 1) or equivalently Z 2 ® SO(2) for N = 2 and 
Z 2 ® SO(N) /SO(N - 1) for N > 2. For XY spins, there exists an identical 
BS as for the STA with first nearest neighbor interaction only. Therefore 
the two systems for N = 2 should belong to the same universality class. 

No numerical studies have been done on this model for the three- 
dimensional case. The two-dimensional case has been studied by Loison 
and SimonPS 

2.3.4. The simple cubic 3\~Ji lattice 

The simple cubic J\—Ji lattice is similar to the stacked J\—J<i square lattices 
shown above, but the second neighbor interactions are also present in the 
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(a) (b) 

Fig. 5. Ground State of the stacked J1—J2 lattice with first (Ji) and second (J2) - 
neighbor antiferromagnetic interactions. J2/J1 > 0.5. (a) and (b) show the two possible 
configurations at non zero temperature. 

xz and yz planes. Similarly, by "order by disorder" the spin configuration 
is collinear, but now there are three ways to choose the direction of parallel 
spins (the z axis in addition to the x and y axes). Therefore the BS is 
Z 3 ® 0(N)/0{N - 1) ee Z 3 <g> SO(N)/SO(N - 1). It is equivalent to the 
BS of the STA with intermediate NNN interaction. 

This model has been studied for Heisenberg spin by Alonso et alp3 and 
Pincttcs and Diep^* 1 . 

2.3.5. J1-J2-J3 lattice 

The addition of a third-neighbor antiferromagnetic interaction to the previ- 
ous model leads to new GS and therefore to a new BS. For some value of the 
interactions a non planar GS state appears!^ In addition a Z 2 symmetry 
is broken from the lattice group and the BS is hence Z 2 <8> 0(N)/0(N — 2) 
between the ordered phase and the disordered phase. For XY spins the BS 
is equivalently Z 2 ®Z 2 &SO(2), and for Heisenberg spins Z 2 ®SO{2). There 
exist also various transitions between various phases, but following the dis- 
cussion at the end of section Ed concerning STA with NNN interaction, 
it is not difficult to find that the transition will be of first order or of the 
ferromagnetic type SO(N - l)/SO(N - 2). 

No numerical studies have been done so far on this model. 

2.3.6. Villain lattice and fully frustrated simple cubic lattice 

The fully frustrated square lattice, 

called Villain lattice, is shown in Fig. EK- This two-dimensional lattice has 
been extensively studied4 107 ll^lf)ff4|llfi|llri|iniD8 119 123. 
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Fig. 6. a. Ground state of the stacked Villain lattice. The double lines are antiferro- 
magnetic interactions, b. The fully frustrated cubic lattice. 

The fully frustrated simple cubic lattice is made from the Villain lat- 
tice in three dimensions (see Fig. |SJd) . The lattice has three ferromagnetic 
interactions and one antiferromagnetic one for each plaquette. In this case 
the GS has an infinite degeneracy for He isenberg spins and is twelve-fold 
degenerate for XY spins. Diep et aP^Eland Alonso et aP^have studied 
this model. 

2.3.7. Face-centered cubic lattice (fee) 

The face-centered cubic lattice (fee) has a degenerate GS. But like the STA 
with NNN interaction, the "most collinear phase" is selected by thermal 
fluctuations. If we consider the smaller unit tetrahedral cell (see Fig.[7Jl, the 
GS consists of two parallel and two antiparallel spin pairs. There are three 
ways to construct. This is equivalent to a Potts Z$ symmetry. Therefore, 
since the GS is collinear, the BS will be Z 3 ®0(N) / 0(N — 1) or, equivalently, 
Z z <g> SO(N) /SO{N - 1). This BS is identical to the STA with intermediate 
NNN interaction. 

This model has been studied by Diep and KawamurgP^ and Alonso et 
alP 

2.3.8. Hexagonal-close-packed lattice (hep) 

The hexagonal-close-packed lattice (hep) has many common features with 
the fee. It is constructed by stacking tetrahedra (see Fig. 01 and an "or- 
der by disorder" mechanism lifts one part of the degeneracy of the GS 
resulting in a collinear GS composed with two parallel and two antiparallel 
spins. Equivalently to the fee case, the BS will be Z 3 <g) 0(N)/0(N — 1) or 
equivalently Z 3 <g> SO{N)/SO{N - 1). This BS is identical to the STA with 
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Fig. 7. a. Ground State of the fee lattice, b. Ground State of the hep lattice. 

intermediate NNN interaction. 

This model has been studied by Dieppe 

2.3.9. Pyrochlores 

The pyrochlore lattice can be made by stacking Kagome lattices along (111) 
direction and is composed of an arrangement of corner-sharing tetrahedra 
(see the chapters by Bramwell et al, and by Gaulin and Gardner, this book). 
In presence of the first and third nearest neighbor interactions Reimera^ 
has proved that the GS is collinear. Similar to the fee and hep cases, there 
are three ways to place the spins on each tetrahedron. The BS should 
therefore be Z 3 ® 0(N)/0(N - 1) or equivalently Z 3 ® SO(N)/SO(N - 1). 
This BS is identical to the STA with intermediate NNN interaction. 
Reimers et alP^ have studied this model. 

2.3.10. Other lattices 

We could stack a Zig-Zag modef^and various phases and phase transitions 
should appear. With an equivalent analysis described at the end of the 
section 12.3.11 concerning STA with NNN interaction, it is not difficult to 
find the nature of the transitions. Still, other lattices can be considered, 
however as shown later, they will have equivalent BS. 



2.3.11. STAR lattices 

In addition one can construct spin systems having an identical breakdown 
of symmetry even if they do not correspond to a real system. 

The first example is derived from the STA. Following the chapter of 
Delamotte et al. in this book ^2 certain modes are irrelevant near the 



February 2, 2008 14:16 



WSPC/Trim Size: 9in x 6in for Review Volume 



loison3 



Frustrated vector spin systems: numerical studies 



13 



critical point. In particular we can construct cells of three spins which are 
always in the ground state 120° configuration. One gets a system of cells in 
interaction but with a rigidity imposed inside the cells (see Fig.|HJ). We call 
this system STAR (Stacked Triangular Antiferromagnetic with Rigidity). 
The BS is then 0(N)/0(N - 2), identical to STA with NN interaction. 

This model has been st udied in t hree dimensions by Dobry and DiepP^ 
and by Loison and Schotte!*SE3I21 

2.3.12. Dihedral lattices Vjv,2 

The second example is derived from the STAR. It is composed of two vector 
spins ei and e2 constrained to be orthogonal to each other at each lattice 
site. The interactions are set to be ferromagnetic and the spin ei (e2) 
interacts only with the other spins ei (e2) at other sites (see the Fig. EJ. 
This model is referred to as the dihedral model Vn, 2- We note that for XY 
spins N = 2 the model can be right-handed or left-handed. 
The Hamiltonian is defined by 



where P = 2. 

At high temperatures the symmetry is O(N) for the first vector and 
0(N — 1) for the second one, since the two vectors must be orthogonal. At 
low temperatures the symmetry is 0(N—1) for the first vector and 0(N— 2) 
for the second one for the same reason. Therefore the BS is 0(N)/0(N — 2), 
identical to the STA with NN interaction. 

This mo del has been studied in three dimensions by Kunz and 
ZumbacfP^I anc | by Loison and Schotte P^^ , and in two dime nsions for 
XY spins by Nightingale, Granato, Lee, and Kosterlitz.^2^22 

2.3.13. Right-handed trihedral lattices 

For Heisenberg spins one can construct another model from the dihedral 
V3 2 with an identical BS. By adding a third vector e3 = ei x e2 to the dihe- 
dral model no degree of freedom is added. Therefore the BS is unchanged 
and identical to the dihedral and to the STA with NN interaction, 



This model has been studied in three dimensions by Loison and Diep. 



p 




(4) 



<ij> k=l 



0(3)/0(l) = 50(3). 
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e(j) 



e(j) 



e,(l) 



(a) 



(b) 



Fig. 8. a. STAR model and dihedral model, b. Trihedral model and V^p model. P = 1 
for J 2 = J-i = 0, P = 2 for J 3 = and Ji = J 2 = 1, and P = 3 for Ji = J 2 = J 3 = 1. 

2.3.14. P-hedral lattices Vn,p 

We can generalize the dihedral Vn,2 model to Vn,p with P orthogonal 
vectors at each site. For P = N the N vectors can be right- or left-handed. 
The BS of symmetry is then 0(N)/0(N - P) with P < N. 

For P — 1 we have the ferromagnetic case, i.e. a collinear GS. For 
P = 2 we have a coplanar GS but no more collinear like in the STA with 
NN interaction. For P = 3 the GS is no more coplanar but restricted to a 
space in three dimensions, and so on. The case N = P = 3 could correspond 
to some experimental systems and spin glasses should also have this kind 
of breakdown of symmetry but in presence of disorder. For P — N the BS 



where Si are iV-component classical vectors of unit length, S qiqj are the 
<7-state Potts spin (q = 2 for Ising spin) with 5 qiqj = when ^ qj, and 
the interaction constant J\ is taken positive (ferromagnetic interaction). 
The sum runs over all nearest neighbors. For N = 2 this model is exactly 
equivalent to the dihedral V^,2 model introduced previously. 

Obviously the BS is Z q ® 0(N)/0(N - 1) or equivalently Z q ® 
SO{N)/SO(N - 1) which is identical to the stacked J x -J 2 model if q = 2, 
and to the STA with NNN interaction, fee, hep if q = 3. 



is O(N)/O(0) = Z 2 ® SO{N). 



This model has been studied by Loison.^ 



2.3.15. Ising and Potts-VN,i model 
We define the following Hamiltonian 




(5) 
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This model has been studied in three dimensi ons by L oison!^ and in 
two dimensions for XY spins by Nightingale et al.^^EZH 



2.3.16. Ising and Potts-VN.2 model 

We can also define a dihedral model with the Hamiltonian JIJ coupled 
with a Potts model in the same way as in J^J. In this case the BS is Z q <g) 
0(N) / 0{N — 2) which is identical to the BS for the stacked J1-J2-J3 model 
if q = 2 and to the STA with a large NNN for q = 3. 
This model has been studied by Loison.E^ 



2.3.17. Landau- Ginzburg model 

The Landau-Ginzburg Hamiltonian is constructed from the dihedral Vn,2 
model (see chapter of Delamotte et al.G2J). We release the constraints on 
the orthogonality of the spins and on the norm unity, replacing them by a 
potential we arrive at the following Hamiltonian 



<ij> i i 

+vY, [(&(*) • M*)) 2 ~ MtfMtf] (6) 

i 

where (f>(i) — (<fi a (i), <fib(i)j is an N + iV-component vector defined on the 
lattice site i, and the summation 53<y> run s over all nearest neighbor 
pairs of the lattice. The first term represents the interactions between the 
sites and the last term the constraint that the spins 4> a (i) and <j>b{i) are 
orthogonal. If v = 0, the Hamiltonian is reduced to a standard 0(2N) 
ferromagnetic problem. This is our third model that we call the chiral (j) 4 
model. For v > the BS is 0(N)/0(N - 2), identical to STA with NN 
interaction. 

This model has been studied numerically by Itakura0 



2.3.18. Cubic term in Hamiltonian 

One can introduce a model where the frustration is not geometrical but 
included on each link resulting in a non collinear GS. Consider the Hamil- 
tonian for two spins: 



H — JiSj.Sj + ^(Si.Sj) 3 



(7) 
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where S, is a N component classical vector of unit length, Ji is a ferromag- 
netic coupling constants J\ < and J3 an antiferromagnetic one J3 > 0. 
This cubic term has the same symmetries as the linear term and for a fer- 
romagnetic system ( J3 < 0) the universality class will be identical!^ For 
1/3 < —J3/J1 < 4/3 the minimum of H occurs when the two spins are 

canted with an angle cos(a) = , 1 On a cubic lattice the GS is pla- 

y/ 3J3/J1 

nar, on a triangular lattice it is in three dimensions (if N > 3), on a fee 
lattice it is in four dimensions (if TV > 4), ■ ■ • This phase has, at finite 
temperature, a transition to a ferromagnetic collinear phase. Therefore the 
BS will be 0{N - 1)/0{N - P) with P = 2 for the cubic lattice (planar 
GS), P = 3 for the stacked triangular, P — 4 for the fee lattice • • ■ The 
case P — 2 has been tested for Heisenberg spins and it gives indeed an 
0(2)/0(l) = 50(2) transition, i.e. a ferromagnetic XY transition.^] The 
case P = 3 has an identical BS as the STA but it has not been tested. 
This model has been studied numerically by Loison.^^ 



2.3.19. Summary 

The number of models that can be constructed by adding interactions is 
unlimited. However, the number of BS can only be finite. As has been 
shown in the previous sections, the symmetry breaking is limited to Z q ® 
0(N)/0(N — P). For the physically relevant cases it is restricted to q = 1 
(Identity), q = 2 (Ising), or q = 3 (Potts), N = 2 (XY spins), N = 3 
(Heisenberg spins), and P = 1, 2 or 3. In total for Heisenberg spins there 
will exist 8 cases (two are identical: (q = 1, P = 3) = (q = 2, P = 2)) and 5 
for XY spins. 

The tables below summarize the probable BS for physical frustrated 
systems with XY or Heisenberg spins. 



3. Phase transitions between two and four dimensions: 
2 < d < 4 

In this section we will concentrate on the nature of the various transi- 
tions mentioned in the previous section. Especially the transition for the 
0(N)/0(N — 2) BS is considered in detail since it appears in numerous 
systems and was extensively debated. Then we will discuss the other BS 
which are less problematic. 
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Table 2. Most probable BS for frustrated XY spin systems with corresponding lat- 
tices or models 



BS 


lattice- model 


SO (2) 


ferromagnetic, small frustration 


Z 2 g> 50(2) 


STA, STAR, y 2 ,2 s= Ising-y 2 ,i, bet, 
Stacked J1—J2, Stacked Villain, 
Fully Frustrated cubic, Stacked Zig-Zag, 
chiral <j> 4 model, (Si.Sj) 3 term 


Z 3 ® SO(2) 


STA+NNN, cubic J1-J2, fee, hep, pyrochlore, Potts- V 2 ,i 


Z 2 8) Z 2 <X> SO(2) 


stacked J1-J2-J3, Ising-V2,2 


Z 3 ® Z 2 <g> SO(2) 


STA+NNN, Potts-V 2 ,2 



Table 3. Most probable BS for frustrated Heisenberg spin systems with corresponding 
lattices or models 



BS 


lattice— model 


SO(3)/SO(2) 


ferromagnetic, small frustration 


SO (3) 


STA, STAR, V 3>2 , right-handed trihedral, bet, 
Stacked Villain, Stacked Zig-Zag, 
Fully Frustrated cubic(?), 
chiral tf> 4 model, (Si.Sj) 3 term 


Z 2 ® SO(3)/SO(2) 


Stacked J1-J2, Ising-V^i 


Z 2 ® 50 (3) 


Stacked J1-J2-J3, ^3,3 =Ising-V3 2 , 
(Si.Sj) 3 term, Fully Frustrated cubic(?) 


Z 2 ®Z 2 ® 50 (3) 


Ising-V3,3 


Z 3 ® SO(3)/SO(2) 


STA+NNN, cubic J1-J2, fee, hep, pyrochlore, Potts- V3,i 


Z- A ® 50 (3) 


STA+NNN, Potts-V 3 , 2 


Z-i <g> Z 2 <g> 50 (3) 


Potts-V 3! 3 



3.1. 0(N)/0(N — 2) breakdown of symmetry 
3.1.1. Fixed points 

bmce many models (STA, STAR, dihedral, chiral </> 4 model, • ■ • , see Ta- 
blcs [2J 01 have an identical BS, they should have equivalent critical behav- 
ior, i.e. they belong to the same universality class. However, the situation is 
more complicated and even with an identical BS two systems could show dif- 
ferent behavior. To understand this fact we have plotted in Fig.Elthe fixed 
points in the critical plan for this model. Since there are two fields (eq. EJ) 
to allow coplanar non collinear GS we will have four possible fixed points. 
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There are two more than the ferromagnetic case which has a collinear GS 
and consequently only one field. 




G H ° G H u G H 

(a) N>N C (b) N = N C (C) N<N C 



Fig. 9. (a), (b), (c) Hamiltonian flow induced by renormalization group transforma- 
tions. The arrows indicate the direction of flow under iterations. 

The fixed points are: 

1. The Gaussian fixed G point at u* = v* = with mean- field critical 
exponents. 

2. The 0(2N) fixed point H at v* = and u* = u H ^ with 0(2N) 
exponents (see Table [JJ. 

3. Two fixed points F + and F_ at location (mj? + ,uf + ) and (mf_ , «f„ ) 
different from zero. These are the fixed points associated with a 
new universality class. 

The existence and stability of the fixed points depend on the number of 
components N: 

a. N > 7V C : four fixed points are present but three are unstable (G, H, 
F_ ) and a stable one F + . Therefore the transition belongs to a new 
universality class different from the standard SO(N)/ SO(N — 1) 
class. If the initial point for the RG flow is to the left of the line 
(G, F-), see Fig.Et, the flow is unstable and the transition will be 
of first order. Therefore two systems with the same Hamiltonian 
and the same breakdown of symmetry could have different critical 
behaviors. 

b. N — N c : the fixed points F_ and F + coalesce to a marginally stable 
fixed point. One would think that the transition is "tricritical" but 
the exponents are different and not given by the tricritical mean- 
field values contrary to common belief. The reason is that there are 
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two non zero quartic coupling constants, see Fig. Et>, in contrast to 
the "standard" tricritical point where the quartic term disappears 
and a sextic term takes over, 
c. N < N c : F_ and F + move into the complex parameter space, 
see Fig. Eh- The absence of stable fixed points is interpreted as a 
signature of a first-order transition. 

There are at least three questions: the value of N c , the location of the initial 
point in the RG flow, and the nature of the transition. 

3.1.2. MCRG and first- order transition 

The most reliable answers to these questions can be found in the article 
of ItakuraP He studied the STA model, dihedral model, and the chiral 4 
model with the Hamiltonian J§J , using Monte Carlo Renormalization Group 
(MCRG). This method is very powerful to give the flow diagram but not 
the critical exponents with a great precision. 

He showed that the STA and dihedral models have an initial point in the 
RG flow under the line GF_ in the Fig. |5| Therefore they should belong to 
the same universality class provided the fixed point F + exists, i.e. N > N c . 
He found that N c is between 3 and 8 in three dimensions which means that 
the real physical systems XY (N — 2) and Heisenberg spins (N = 3) have 
a first-order transition. 

In addition he did a standard canonical Monte Carlo (MC) simulation 
for XY spins on the STA lattice for very large sizes of 96 3 and 126 3 . He 
found a first-order transition in agreement with MCRG study, co ntrary to 
smaller size systems which seem to have a second-order transition!™^ For 
Heisenberg spins he could not find a first-order transition but using the 
results of the MCRG he concludes that the first-order transition could only 
be seen for a size larger than 800 3 which is not accessible for actual com- 
puter resources. For the dihedral model with a canonical MC he found a 
clear first-order transition for large sizes for Heisenberg spins. For XY spins 
Loison and Schotte have already shown that the transition is of first order. 

in 1989 was the first to find a first-order transition in helimagnetics 
bet lattice with XY spins. This result was considered not conclusive be- 
cause of problems of periodic boundary conditions in numerical simulation. 
However, as noted above, this is not relevant for XY spins (contrary to 
Heisenberg spins) and whence this conclusion is indeed correct for this BS. 

We have now to address the problem why phase transitions appear con- 
tinuous for small sizes but show the true first-order nature only for larger 
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sizes. This phenomenon is not restricted to frustrated spin systems such 
as STA but appears also in the weakly first-or der trans ition of the two- 
dimensional Potts model with q — 5 components ) 42 | 43 J44| M ore generally it 
appears when two fixed points collapse and disappears following one vari- 
able (like i 7 ! and F + in Fig. EJ) . See later for a comparison with the Potts 
model. 

3.1.3. Complex fixed point or minimum in the flow 

The change from continuous to discontinuous transition can be understood 
using the concept of "c omp lex fixed point" or "minimum in the flow " first 
introduced by ZumbacrP^and then developed by Loison and Schotte.^ESI 
In Fig. for N < N c the fixed points F_ and F + have collapsed and no 
solution exists in real parameter space. Instead there exists an imaginary 
solution plotted in Fig. HUb . These solutions should have an influence on 
the flow in the real plane as shown in Fig.llOb. Zumbach showed that there 
exists a basin of attraction due to the complex fixed points where the RG 
flow is very slow. 



real space 


ttk H 


(u,v) Q 


F 




complex axe 


(a) 


N<N C 




Fig. 10. a) For N < N c the fixed points _F+ and F— become complex, b) Hypothesis on 
the Hamiltonian flow induced by renormalization-group transformations for TV < N c (see 
Fig-IEt)- The arrows give the direction of flow. The two circles correspond to a low velocity 
region. The inner one corresponds to a minimum and hence an "almost" second-order 
transition. The outer circle corresponds to a slow crossover where the critical exponents 
will vary as a function of the system size. 



If the RG flow of a system goes through this region, the number of RG 
iterations L — > L/b to get outside will be large. Therefore the size L of the 
numerical system considered, should be made large enough (seeEJ. If the 
size is "too small" , the flow is trapped inside of this low velocity region 
and the transition seems to be continuous. We immediately see that the 
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size L for which the true first-order transition is visible, that is when the 
flow is outside of the basin of attraction, will depend on the starting point. 
If it is located outside of the domain of low velocity and if the flow does 
not approach this region (point A in Fig. IKljl , then the first-order nature of 
transition can be seen already for small sizes. In the other case, if the initial 
position is such that the flow has to go through the entire region of low 
velocity (point B in Fig. ll(J|) . then the necessary size to reach the true first- 
order region will be very large. This explains why STA, STAR, dihedral and 
right-handed trihedral (for N — 3) models having an identical breakdown 
of symmetry behave differently. Similarly a single system can show different 
behaviors for different sizes. For the Hciscnbcrg case, for example, the right- 
handed trihedral model shows a strong first-order behavior even for small 
sizes. The dihedral model with the same BS for N — 3 has a weaker first- 
order behavior visible only for sizes bigger than 80 3 spins while the STA or 
bet helimagnets show a second-order transition for similar sizes. 

Furthermore the size of this region of low velocity will vary as a function 
of the distance between the complex fixed point F + and the real plane shown 
in Fig. HOfa. It has been shown that N c , where the two fixed points F_ and 
F + collapse and become complex, is bigger than 3 and the distance of F + 
to the real plane will be larger for N = 2 than for N = 3. Therefore the 
size of the domain of low velocity will be bigger for N — 3 than for N = 2, 
and the true first-order behavior will be seen for smaller sizes in the XY 
case (N — 2) than in the Heisenberg case (N — 3). Indeed the first-order 
transition in the dihedral, STAR and STA can be seen from the size 12, 
18 and 96, respectively, for XY spin systems (see Table 0J). For Hciscnbcrg 
spins this has been seen for the dihedral model with a size L = 80 (see 
TableEJ) while the estimate size for the STA should be larger than 8OO0We 
note that a first-order transition has been found in quasi-one-dimensional 
STApS the initial point in the flow diagram is outside the domain of low 
velocity. 

For the XY case Delamotte et al@^ using a non perturbative RG ap- 
proach, found that there is no minimum in the flow but always a low ve- 
locity region (outer circle of Fig. II (lb). Following the initial point in the 
RG flow the systems will show different sets of exponents, i.e. the system 
is in a crossover region. Numerical simulations tend to support this in- 
terpretation. Indeed for XY spins the STAR and the STA for small sizes 
(L < 40) display a second-order transition with different critical exponents 
(see Tables EH5)l and for large sizes the STA shows a first-order transition 
as discussed above0 
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In conclusion the transition for Heisenberg spins is of first order but an 
"almost second-order transition" could exist for a range of "finite" sizes. 
This "almost second-order transition" will have a set of exponents different 
from the ferromagnetic ones because the breakdown of symmetry is differ- 
ent: 0(N)/0(N — 2) in comparison to 0{N)/0(N — 1) for ferromagnetic 
systems. This "new" universality class has been called chiral class For the 
XY case there is a crossover region with a slow velocity and exponents will 
vary between those of the chiral class and of weak first-order transitions 
{v = 1/3 - see Table QJ. The critical exponents of this chiral class are 
given in Tables EE for N varying from 2 to 6 for the STA, STAR, and 
dihedral models. With numerical simulations, exponents j/u (3/v and v 
are usually calculated using finite size scaling while other exponents are 
calculated using the scaling relation dv = 2 — a and j/u = 2 — r\ (seeJBJ. 

There is another indication given by Loison and Schott J21 which sug- 
gests that the phase transition is indeed due of a "complex fixed point'. 
Using the scaling formula 7/1/ = 2 — 77, one gets a negative r\ for N = 2 and 
N = 3 for the STA, STAR or dihedral model (see Tables [TO which con- 
tradicts the theorem that 77 must be always posi tive! V ' ' ' Negative 77 can 
also be calculated using experimental results!^ Therefore the fixed point 
cannot be real. Using this indication we can conclude that the continuous 
transition found in fully frustrated XY lattice- is indeed of first order. 

Table 4. Critical exponents associated to the SO (2) symmetry by Monte Carlo for 
XY spins (N = 2) and a BS Z 2 ®SO(2). < a ) Calculated by f/v = 2 -t}. The first result 
comes from a study at high and low temperatures and uses the FSS. The seconc0uses 
the Binder parameter to find T c and uses the FSS, the thircJZIuses the maxima in FSS 
region. 



system 


Ref. 


Lmax 


a 


P 


7 


V 


V 


STA 


m 


60 


0.34(6) 


0.253(10) 


1.13(5) 


0.54(2) 


-0.09(8) (a) 


STA 


n 


33 


0.46(10) 


0.24(2) 


1.03(4) 


0.50(1) 


-0.06(4)( a ) 


STA 


□ 


36 


0.43(10) 






0.48(2) 




STA 


HI 


126 


first order 


STA 


m 


35 


first order 


bet 


m 


24 


first order 


STAR 


El 


36 


first order 


V 2 ,2 


M 


36 


first order 
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Table 5. Critical exponents associated to the Z2 symmetry (chirality n) by Monte 
Carlo for XY smns (N = 2) and a BS Z 2 ® SO (2). O) Calculated by j/u = 2 - J7. 
The first result "I comes from a study at high and low temperature and uses of FSS. 
The second Ifil uses the Binder parameter to find T c and uses the FSS. 



system 


Ref. 


Lmax 


ct 




Ik 






STA 


HI 


60 


0.34(6) 


0.55(4) 


0.72(8) 


0.60(3) 


0.80(19)0) 


STA 


M 


33 


0.46(10) 


0.38(1) 


0.90(2) 


0.55(1) 


0.28(3)0) 


STA 


El 


126 


first order 


STA 


M 


35 


first order 


bet 


19 


24 


first order 


STAR 


IMI 


36 


first order 


V 2 ,2 


IMI 


36 


first order 



Table 6. Critical exponents by Monte Carlo for Heisenberg spins (N = 3) and a BS SO(3). 
Calculated by O) j/u = 2 - rj , 0) du = 2 - a , O) 2 /3/u = d - 2 + r) . 



system 


Ref. 


Lmax 


a 


/3 


7 


V 


V 


STA 


m 


60 


0.240(80) 


0.300(20) 


1.170(70) 


0.590(20) 


+0.020(180)0) 


STA 


m 


36 


0.242(24)0) 


0.285(11) 


1.185(3) 


0.586(8) 


-0.033(19)0) 


STA 


IT2l 


48 


0.245(27) < 6 ' 


0.289(15) 


1.176(26) 


0.585(9) 


-0.011(14)0) 


STA 


rrn 


36 


0.230(30) (6) 


0.280(15) 




0.590(10) 


0.000(40)0) 


bet 


m 


42 


0.287(30)( fc ) 


0.247(10) 


1.217(32) 


0.571(10) 


-0.131(18)0) 


STAR 


m 


42 


0.488(30) 0) 


0.221(9) 


1.074(29) 


0.504(10) 


-0.131(13)0) 


V 3 ,2 


1551 


40 


0.479(24)0) 


0.193(4) 


1.136(23) 


0.507(8) 


-0.240(10)0) 


V 3 ,2 


E 


80 


first order 


Table 7. Critical exponents by Monte Carlo for spins with four components (N = 4) and a 
BS 0(4)/0(2) = SO(4)/SO(2). Calculated by < a ) j/u = 2 - n , < 6 ) du = 2 - a . 


system 


Ref. 


| Lmax 


I a 


1 P 


1 7 


] v 


V 


STAR 


I5B1 


42 


0.287(27)( 6 ) 


0.291(11) 


1.133(28) 


0.571(9) 


+0.015(18)0) 


V 4 ,2 


m 


40 


0.278(30)( 6 ) 


0.290(12) 


1.142(34) 


0.574(10) 


+0.011(25)0) 


Table 8. Critical exponents by Monte Carlo for spins with six components (N = 6) and a 
BS 0(6)/0(4) = SO(6)/SO(4). Calculated by < a ) y/v = 2 - n , ( b > du = 2 - a . 


system 


Ref. 


Lmax 


a 


1 


I 7 


1 V 


I v 


STA 




36 


-0. 100(33) W 


0.359(14) 


1.383(36) 


0.700(11) 


+0.025(20)0) 
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3.1.4. Experiment 



A resembling situation occurs in the analysis of experiments. There the 
correlation length £ plays the role of the system size in numerical simu- 
lations. For a second-order transition one has £ ~ (T — T c )~ v . Therefore 
we should observe a crossover between the region of low velocity ( "almost 
second-order transition" in Zumbach's words) to the true first-order behav- 
ior for temperatures "closer" to the critical temperature. Unfortunately the 
situation is even more complicated in experimental systems with the om- 
nipresence of planar or axial anisotropies and the one- or two-dimensional 
characters of the compounds. Then a succession of crossovers (sec [6] for 
more details about crossovers) from 2d to 3d and from Heisenberg to Ising 
or XY behavior could lead to difficulties in the interpretation. Furthermore 
other (small) interactions could also dominate the behavior near the critical 
temperature and change the universality class. 

Experiments for XY spins: All experimental results can be found 
in Ref. [34]. Several A XB^ compounds have the STA structur e. The ex- 
periments on CsMnBr 3 | 4» l ^) | 51 | 52 | 5:i | ^4 | 55 | .^ | 57 | 58 | ^ | RbMnBr Jilll and 
CsVBrJSS 

give critical exponents compatible with those of MC simulation 
on STA and a second-order transition. We can interpret this result by the 
fact that the systems are under the influence of a complex fixed point, and 
t oc T — T c is too sma ll to observe a first-order transition. 

The case CsCuClJ^l is different since the authors observe a crossover 
from a second-order region with exponents compatible with MC results 
on STA for 10~ 3 < t < 5.1CP 2 to a region of first-order transition for 
5.10 -5 < t < 5.1CP 3 . For t < to w 10~ 3 one seems to observe the true 
first-order region which corroborates the scenario introduced previously. 

Three kinds of helimagnctic structure have been studied: Holmium, Dys- 
prosium and Terbium. Essentially three types of results exist: those com- 
patible with MC ones of the STA, those with a large (3 incompatible with 
STA and those showing a weak first-order transition. 

The results co mpatible with those of MC on STA for Ho ES1MD 
and To ' ^ l au l ° 1 l sz l can be interpreted as before: the sys- 
tems a re und er the influence of F + . The first-order transition for Ho«^l 
and DjE^ZH i s d ue to the fact that the measurements were done in the 
first-order region near the transition temperature. The values of the ex- 
ponent (3 ~ 0.39 in the case of Hcl 67 ! 68 ! 69 ! 7 ^ and DyEHESE3Z3 are not 
compatible with those found by MC ((3 ~ 0.25). This fact can be explained 
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by the presence of a second length scale in the critical fluctuations near T c 
related to random strain fields which are localized at or near the sample 
surface!^ Thus the critical exponent j3 measured depends on this second 
length. 

Experiments for Heisenberg spins: All experimental results can 
be found in Rcf. [34]. As explained previously, before the first-order 
region is reached, the crossover from Heisenberg to Ising or XY be- 
havior prevents a first-order transition of Heisenberg type. Neverthe- 
less the second-order transition can be studied for VCla P ^ VBrJ^l 
Cu(HCOO) 2 2CO(ND 2 )22D 2 dHH an d Fe[S 2 CN(C 2 H 5 ) 2 ] 2 ClPlFor the last 
two examples the observed exponents might be influenced by the crossover 
from 2d to 3d Ising behavior. The experimental results agree quite well with 
the MC simulations. 

3.1.5. Value of N c 

Tables 14181 allow us to get an idea of the value of N c . The MCRG gives 
an estimate 3 < N c < 8 but the values of ij can give a better estimate 
following Ref. [34]. As can be seen for the XY and Heisenberg systems, 
negative values of t) ap pear fo r the STA, STAR or dihedral model. But r\ 
must always be positive { 46 | 47 | ^his j s c [ ue ^ th c usc of the scaling relation 
7/V = 2 — r\. Indeed ZumbacfP^ has shown that for an "almost second- 
order transition" , i.e. when the solution becomes complex, this relation has 
to be modified to 7/V = 2 — ry + c, c being a constant different from zero. 
We can use this relation as a criterion for real or complex fixed points. In 
three dimensions rj is usually small and is almost independent of N for the 
ferromagnetic case, that is ~ 0.03 (see Table^). Our hypothesis is that it is 
also true for the frustrated case. Indeed for N = 6 we found r\ ~ 0.03 (see 
Table|HJ). Accepting this value c becomes zero around N c ~ 4.5. Obviously 
if a bigger value of rj is chosen, N c will increase. 

3.1.6. Phase diagram (N,d) 

In Fig. 1111 we have plotted a phase diagram where the abscissa is the spin 
dimension N and the ordinate the space dimension d. 

There is a line which divides a region of first-order transition from a 
region of second-order transition. At, and near, d = 4 we can use the RG d — 
4 — e expansion to find N c — 21.8 — 23.4 e, and for d — 3 we have seen that 
iV c ~ 4.5. Furthermore there is a region near this line where the transition is 
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Fig. 11. The critical curve N c (d) separating the first-order region for small N and large 
d from the second-order region for large TV and low d. The squares represent the physical 
systems of interest (N = 2, 3 in three dimensions). 



"almost second order" , i.e. under the influence of the complex fixed point for 
"small" sizes. The system will display a second-order transition with stable 
critical exponents. Besides, there is a slow crossover region where critical 
exponents vary with the system size. The sizes of these regions depend on 
the initial point in the flow diagram and therefore on the model. 

Near two dimensions and for Hciscnbcrg spins the transition is of the 
type SO(4)/SO(3), i.e. ferromagnetic for spins with four components (see 
the following section and the chapter of Delamotte et alPZlfor more details). 



3.1.7. Renormalization- Group expansions 

We present here a review of the RG expansions ( "non perturbative" , d = 
4 — e, in fixed dimension d = 3, and d = 2 + e). For more details see El and 
the chapter of this book of Delamotte et al!"*l 

We discuss now the 4 — e expansion. Using a continuous limit of the Lan- 
dau Ginzburg Hamiltonian it involves a perturbative extension (u,v) 
around the Gaussian solution with the dimension d and the number of spin 
components N as parameters. To get the result for d = 3 one puts e = 1. 
Of course the series are at bes t asymptotic and must be resummed, even 
for the ferromagnetic case !^ The extension until e 2 in Ref. [88] has been 
"resummed" (only three terms) and the value of N c (d = 3) ~ 3.39 is com- 
patible with the numerical result _/v™ ume ™ cai (rf = 3) ~ 4.5. The calculations 
have been extended to e 4 in Ref. [89] but unfortunately no resummation 
has been done on this series yet. 

The expansion for fixed dimension d = 3 is very similar to the d — 4 — e 
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one. The series and the resummation will differ slightly. At three loops the 
result is N c ~ 3.9 J^which is again not far from N? umerical (d = 3) ~ 4.5. A 
very interesting picture emerges from the six-loop calculations! ^ ' First, 
for N ~ 5, the resummation does not converge. Second, for TV = 3, the 
resummation gives rise to a very large variation of the critical exponents 
following the flow chosen to arriveat a single critical exponent at the critical 
point (see Fig. 5 of the referencd^J) which is a really "new" phenomenon. 
Third, for N = 2 a second-order transition is predicted in contradiction 
to Itakura's result JS1 and the numerical simulations (see Tables PH5|l . Man- 
ifestly we could conclude that the resummation chosen does not work for 
this series. To understand this we have to stress that the series has many 
more terms because of the double expansions in u and v compared to a 
single one for the ferromagnetic case. Whence it is much more difficult to 
find a "good" resummation scheme. 

On the other hand Delamotte et alP^ found that the transition must 
be of SO(A)/SO(3) type, i.e. equivalent to a ferromagnetic 0(4) transition, 
near two dimension for Heisenberg spins. This study was done using the 
d = 2 + e expansion or equivalcntly the non linear a model and the result is 
valid only nea r two dimensions. This prediction was checked by Southern 
and Young® by MC in two dimensions. But we can give argumentJ^ 
using l/W-expansion and the value of the critical exponents to rule out 
this kind of transition in three dimensions. These arguments hold as long 
as the relevant operators are identical in l/N expansion and in the non 
linear a model. Indeed for the l/N expansion, and equivalently for the 
4 — e expansion and in fixed dimension, we keep in the Hamiltonian only 
the terms which are renormalizable near d = 4 or infinite N, and discard the 
others. Nevertheless it seems that some non renormalizable terms (seeEJ 
become relevant and important for low N and d and, whatever the number 
of loops we cannot find the correct behavior using these expansions. 

Because of the problems encountered by the standard expansions, a "non 
perturbative" approach could be very useful. This was introduced in this 
model by Tissier et al.ESI We quote "non perturbative" because it is not an 
expansion resembling the other methods where the extension parameters 
(u, v, e) are usually not small. In this non perturbative method even if we 
introduce only a few terms in the action, results will be very good and no 
resummation is necessary. We notice that the simplest action was studied 
by ZumbachE£l which allows him to introduce the "almost second-order" 
transition. Adding more terms, Tissier at al. found that they are able to 
retrieve all the previous expansions with additional information. Their value 
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of iV c ~ 5.1 is comparable to N™ umencal ^ 4.5. i n three dimensions they 
found critical exponents very close to those calculated by MC for N = 6 (see 
Table|BJl, a minimum in the flow for N = 3, and a slow crossover for N = 2 as 
explained in the previous section. They found also that the critical behavior 
is indeed that of a ferromagnetic 0(4) transition near two dimensions for 
Heisenberg spins. Some non-renormalizable operators excluded in the 4 — e 
and 1 /N expansions are included in this "non perturbative" method. These 
operators are always relevant between two and four dimensions whatever 
N is, but have an influence on the values of the critical exponents only for 
small N and near two dimensions. This explains the discrepancy between 
the d — 2 4- e and d = 4 — e expansions. 



3.1.8. Short historical review 

In this section we give a short historical review of studies on frustrated 
systems. Indeed the history was not straight if we look back at the last 25 
years. The de velopment began by a RG 4 — e expansion by Jones, Love 
and MoorJSI in 1978. It is only in 1984 and in the following years that 
Kawamurg0 started with the first numerical simulations and, in combina- 
tion with the results of several experiments, proposed a "new" universality 
class. 

Then Azaria, Delamotte et ai^l found that the transition near two di- 
mensions for Heisenberg spins should belong to the 0(4) ferromagnetic 
class. Several groups (see Tables l4T5)l have done numerical simulations and 
experiments which favored either the new universality class, a first-order 
transition, an 0(4) transition a nd eve n a mean-field tricritical transition. 
During those years Sokolov et al.^3^ have extended the RG expansion to 
three loops and found a first-order behavior for XY and Heisenberg spins. 

Surely one of the most important articl es to find the key to under- 
stand the physics was written by ZumbachpSI using a non perturbative 
approach. He predicts an "almost s econd- order transition" for XY and 
Heisenberg spins. Loison and SchottepEH using this concept, were able to 
get a clear picture for both the numerical and experimental studies. Then 
Tissier, Delamotte and Mouhanna by extending the work of Zumbach, 
were able to understand the whole phase diagram for the dimension be- 
tween two and four. To terminate Itakura,^ confirmed the picture given by 
Zumbach&Loison&Schotte, which we think is the definitive answer from a 
numerical point of view. 
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3.1.9. Relations with the Potts model 

We would like now to stress the similarities between the Potts model and 
the frustrated case just studiedESIln Fig. ^| we have plotted the RG flow 
diagra rrP as a function of the first- and second-neighbor ferromagnetic 
interactions (Ji and J2) and the chemical potential A (corresponding to 
the site vacancy). 




(a) t<t c (b) q=q c (c) q>q c 

Fig. 12. Hypothesis of the renormalization flow in presence of second-neighbor inter- 
action. For q < q c the initial point of the flow is irrelevant to the critical behavior. For 
q > q c the behavior for accessible lattice size will be different for the curve 1 (going 
through the basin of attraction of the fixed points F+) and the curve 2 (going directly 
to the first-order fixed point D). 

F + is the standard stable ferromagnetic fixed point and F_ is an unsta- 
ble fixed point. D is a first-order fixed point. O is the initial point of the 
system. This figure resembles the one of the frustrated case (Fig. [5J). For 
q < q c the flow goes to the fixed point F + and has a second-order transition. 
At q = q c = 4 (in two dimensions) the two points F + and F_ collapse and a 
second-order phase transition with logarithm correction results. For q > q c 
there is no longer a solution and the flow goes directly to the fixed point D, 
i.e. the system has a first-order transition. Our hypothesis is that even if the 
solution of F+ is complex, it has an influence on the real space and there is 
a region of low velocity in the flow diagram. It can be a true minimum for 
q > q c with an "almost" second-order transition. For q bigger, no minimum 
should exist but, following the sizes studied, a variation of the critical expo- 
nents could occur. Indeed for q — 5 a set of critical exponents was found™" 
corresponding to our hypothesis of "almost second-order" transition. The 
cases q = 6 and q = 7 have not been studied, but possibly one could observe 
a variation of critical exponents corresponding to a slow crossover to the 
first-order point. 
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We remark that depending of q we can observe clearly if the transition 
is of first order or second order. In Fig. ^| we have plotted the flow (2) in 
presence of the second-neighbor interaction. We observe that this flow is 
less influenced by the low velocity region than the flow (1) if q > q c and it 
reaches F + if q < q c . Adding the third- nearest-neighbor interaction, we are 
able to find numerically that q c = 4 for the crossover from the first-order to 
second-order transitions. Contrary to the frustrated cas e we know roughly 
the coordinate of the fixed point F + (J 2 /Ji ~ 0.3^3 an d A ^ O^By 
consequence we know which kind of interaction must be added. 



3.2. O(N) /0(N — P) breakdown of symmetry for d — 3 

In mean-field theory, for N > P the model shows a usual second-order 
type, but for N — P the transition is a special onel-^The BS in this case is 
Z2®SO(N) and the coupling between the two symmetries leads to a special 
behavior which, for N = P = 2 in two dimensions, has been extensively 
debated (see later in this chapter). 

These generalized chiral mode ls have b een studied by applying the RG 
technique (d = 4 — e expansion) pv l ytf l vV I The picture is very similar for 
all N > P > 2. At the lowest order in e, there are up to four fixed points, 
depending on the values of N and P. Amongst them arc the trivial Gaussian 
fixed point and the standard isotropic 0(NP) Heisenberg fixed point. These 
two fixed points are unstable. In addition, a pair of new fixed points, one 
stable and the other unstable, appear when N > N c (d) with 



N c (d) = 5P + 2 + 2v/6(P + 2)(P- 1) 
25P 2 + 22P - 32 

5P+2 



2 V /6(P + 2)(P-1) 



(8) 



For P = 2 we find the standard result N c = 21.8 — 23.4 e. On the other 
hand, for P = 3 we obtain N c = 32.5 — 33.7 e and for P = 4 we obtain 
N c — 42.8 — 43.9 e. A "tricritical" line exists which separates a second-order 
region for low d and large from a first-order region for large d and small 
N. Applying e = 1 (d = 3), we obtained that N c (d = 3) < for all P. 
For P — 2 we know that this result does not hold and equivalently it does 
not apply for P > 3. Loison 03 has done some simulations for N = 3 and 
N = 4 with P = N and P = N - 1 for the Stiefel V N , P model. He showed 
that the transition is clearly of first order using numerical simulations. He 
generalized this result for all N . We remark that the fully frustrated cubic 
lattice (generalized Villain lattice) which could have an identical BS has also 



February 2, 2008 14:16 



WSPC/Trim Size: 9in x 6in for Review Volume 



loison3 



Frustrated vector spin systems: numerical studies 31 

a first-order transition for Heisenberg spinsP^In conclusion the transition 
in systems with a O(N) / 0(N — P) the breakdown of symmetry is of first 
order for P = N and P = N — 1 and in particular for the physical case 
N = P = 3, i.e. a non planar GS for Heisenberg spins. 

3.3. Z 2 <g) SO(N)/SO(N - 1) breakdown of symmetry for 
d — 3 

For N > 2, this BS corresponds to a collinear ground state associated to a 
BS of the lattice (see Table I""!""*) . For example the stacked Villain lattices 
with first-neighbor ferromagnetic interaction J\ and a second-neighbor anti- 
ferromagnetic interaction J2 can have this BS. For J2/J1 < —0.5 the compe- 
tition between the two interactions leads to two possible ground states!^ 
For N = 2 (XY spins) the stacked Villain lattices and the stacked J\- 
J2 model have the same BS as the STA and should display an analogous 
behavior ("almost second order" for small sizes following by a first-order 
transition for larger sizes). 

We have simulated by MC technique the stacked Villain lattices and 
the stacked J1—J2 lattice which have this kind of BS. We found that the 
transition seems continuous^"! for small lattices for XY and Heisenberg 
spins. For XY spins this is in agreement with the numerical simulation of 
STA which has an identical BS. We have also studied the Ising— Vjy,i model 
which has the same BS. For N — 2 this model is the equivalent to the V-^2 
model which has a strong first-order behavior as explained previously. For 
N = 3 and N = 4 we found also strong first-order transitions. Therefore it 
seems that the transition Z2 <£> SO(N)/ SO(N — 1) is of first order for any 
N even if it could exist an "almost second-order transition" for small sizes. 

3.4. Z 3 (g) SO(N)/SO(N - 1) breakdown of symmetry for 
d — 3 

This BS appears with a collinear ground state associated to a three-state 
Potts symmetry due to breakdown of lattice symmetries. This corresponds 
to various lattices where the unit cell is composed of four spins (two par- 
allel spins and two antiparallel spins) like the STA with intermediate NNN 
interaction (0.125 < J2/J1 < 1), the fee, the hep, the cubic J1-J2 model or 
the pyrochlore (see Table 12*131) . 

As seen previously, even the coupling of an Ising symmetry (Z2) with 
SO(N)/SO(N — 1) can give a first-order transition. Whence it is expected 
that the transition with a three-state Potts symmetry (Z 3 ) gives also a 
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first-order transition which is even stronger if we consider that the Potts 
symmetry alone has a first-order transition in d = 3 for q = 3. Indeed a 
stron g first- order transition was seen in the fccp^EHI h cpj^j cubic J\-J% 
mode l 25 ! 26 ! pyrochlorePSI and STA with NNN interactiorCEIl for XY and 
Heisenberg spins. 

We have done some simulations for the Potts- Vn,i model with the 
Hamiltonian (JSJ) with q > 3, i.e. the g-state Potts model coupled to an 
iV-component vector. We found a strong first-order transition for any iV 
and q. 

In conclusion the transition for the Z q ®SO(N)/SO(N -1) BS in d = 3 
is of first order for any N > 2 and q > 2. 

3.5. Z q ig) 0(N)/0(N - 2) and other breakdown of 
symmetry in d — 3 

This BS appears with a planar ground state associated to a three-state 
Potts symmetry. One example is the STA with NNN interaction J2I J\ > 1 
(see Table EP. 

To simulate this BS we have used a Potts-V/v^ model for any q. We 
found the behavior of a strong first-order transition whatever N and q are. 
It is not surprising if we consider that the "less frustrated" model Potts- 
Vn,i is already of first order. 

We have also simulated the Z q ® 0(N)/0(N - 3) for N = 3 which 
corresponds to the BS Z q ® Z 2 ® SO{N) using a Potts-V/v,p model. We 
found again a strong first-order transition. Similar results are obtained for 
Ising-Potts-Vjv,! and Ising-Potts-V/v,2 models. 

To summarize, it seems that the coupling between a Potts or Ising sym- 
metry and a continuous symmetry of vector spins give always rise to a 
first-order transition. 

4. Conclusion 

We have studied breakdowns of symmetry of three-dimensional lattices for 
XY (N — 2) and Heisenberg (N = 3) spins. The general form of the break- 
down of symmetry is Siattice ® 0(N)/0(N - P) . Si aU ice is usually a Potts 
symmetry Z q (Z\ corresponds to the identity, Z2 to an Ising symmetry, Z% 
to a three-state Potts symmetry, . . . ). P runs from 1 to 3. 

If the frustration is small the GS can be collinear (P = 1) with- 
out breaking any symmetry of the lattice Siattice = 1, the BS will be 
SO(N)/ SO(N — 1) and the transition will be of second order like in ferro- 
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magnetic systems. 

If the frustration is strong enough, we can have a non collinear GS with 
the usual planar configuration (P = 2) like in the STA or helimagnets 
and Siattice = 1- Depending on the model and for not "too big" sizes, the 
transition could appear continuous belonging to a new universality class. 
For an infinite size the system will have a first-order transition. 

Some more complicated systems could also have a non planar GS, and 
for Heisenberg spins it corresponds to P = 3 with usually Siattice — 1< 
The transition will be of first order with a possible "almost second-order" 
behavior for small system sizes similar to the P = 2 case. 

Some even more complicated systems could have a non collinear GS 
(P = 2 or P = 3) with Siattice = Z q and q = 3. For example the STA with 
large NNN interaction has a strong first-order transition. 

On the other hand some frustrated systems have degenerate GS but by 
"order by disorder" a collinear GS (P — 1) is selected. The lattice symmetry 
could be broken giving an additional BS. Siattice = Z-x (Ising) for the stacked 
Villain model, or Siattice = Z% in fee, hep, STA with intermediate NNN 
interaction, and pyrochlore. For the Z% symmetry the transition is strongly 
of first order even for small sizes. For the Zi symmetry the transition is also 
of first order for large sizes but looks continuous for small sizes. Therefore 
we cannot exclude an "almost second-order" behavior belonging to a new 
universality class. 

In conclusion frustrated systems have a first-order transition for XY 
and Heisenberg spins in three dimensions even if for "small" sizes the sys- 
tems would show a second-order transition. "Small" could mean sizes of 
thousands of lattice constants and the first-order transition will then not 
be observable with the actual computer resources. 

5. O(N) frustrated vector spins in d — 2 
5.1. Introduction 

The situations in two-dimensional systems are different from the three di- 
mensions. The Mermin-Wagner theorerrP^ asserts that no magnetiza- 
tion appears at n on-zero temperature . However, a transition due to the 
topological defect can appear. The binding-unbinding transition 
of vortex-antivortex pairs for XY spin systems is a classical example. We 
note that their e xact role for three-dimensional phase transitions is not 
clear UHa fll U4 | 1 U5 | p or non co rj mear Qg ; the topological properties of the 
system differ from those of collinear GS and new types of transition could 
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appear. In the following we will briefly review the ferromagnetic case, then 
the frustrated XY case, applicable to Josephson-j unction arrays in mag- 
netic fields, and finally the frustrated Heisenberg cases. 

5.2. Non frustrated XY spin systems 

For a non frustrated XY spin system the order parameter is 5*0(2). The 
topological defects of this group is only the point defect classified by the 
homotopy gr oup TLt(S O(2)) = Z, Z being the topological qua ntum number 
of the defect .IMLllSl A Kosterlitz-Thouless (KT) transition EH exists at 
the critical temperature T c ~ 0.9, driven by the unbounding of vortex- 
antivortex pairs. 

This transition has some special features: for T < T c the correlation 
length (£) is infinite, while for T > T c it has an exponential decreasing 
contrary to the power-law behavior in the standard transition. 

Several numerical methods exist to calculate the critical temperature 
and the critical exponents (seeEJ). 

5.3. Frustrated XY spin systems: Z?, (g) 5*0(2) 

A frustrated XY spin system can have many order parame- 
ters as previously seen (see Table |2Jl. The most studied case is 
Zi % S O (2) which corr esponds to a non collinear GS: triangular 

Villain . lattice or two-d imensional fully frus- 
trated latticJ 1 '™ 1 ^ lfll U\i lfi|l 1 71113111 

and related Zig-Zag 

models.^^201 i n combination to the topological defect point Z, there is a 
line of Z2 corresponding to an Ising transition. Therefore they will have sim- 
ilar properties as the J\—Ji model, the T/2,2 model or any model where the 
Z2 symmetry comes from the lattice. We note that different models can have 
an identical group and therefore the same topological defects. The following 
model s sho uld have analogous properties : the Vil lain model on the Vi llain 
latticePUl The Ising- VVi = V 2 , 2 modelPSffill the 19- v ertex model 
the ID quantum spin s! _ the Coulomb gas repres enta tion! LdSU^JlLdillLdja 
the XY - XY modelpnUIIHor the RSOS modelPSl 

This kind of models was widely studied because it is believed to 
correspond to experimental systems such as Joseph son- Junction arrays 
of weakly coupled superconducting i 

siandjnmsn 

or films of Helium 

3 He paMOTH] 

The question is to understand the coupling between the vortex- 
antivortex and the Ising symmetry. Contrary to the three-dimensional case 
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discussed previously there are no clear accepted answers to these questions. 
We will therefore only present the different results in the high temperature 
region and in the FSS, giving some clues to understand the physics of theses 
systems. 

What are the possible scenarios? 

The KT transition associated to the topological defects appears at 
a temperature lower than the one associated to the Ising transition 
even if they can be very close. The two transitions should have 
standard behaviors (v KT = 0.5, v lsing = 1, • • • ). 
Furthermore Olsson has proposed that the Ising behavior could be 
observed only when L 3> £ KT '. Since £ KT is infinite for T < T^ T 
and decreases exponentially for T > T c , it will correspond to a high 
temperature region with T > Tf T , but also T > T c /sm9 when 
the two transitions appear in close neighborhood. In particular the 
FSS region will have a non standard Ising behavior for the sizes 
accessible with actual computer resources. Only for a very large 
lattice the standard Ising behavior will be visible. We observe the 
similarity of this hypothesis with the one described previously in 
three dimensions and we could call the transition in the FSS an 
"almost new Ising-KT behavior" . 

2 rpKT ^ rplsing 

The two transitions appear at the same temperature and display 
different behaviors compared to the standard one. We remark that 
for high temperatures T S> T c , i.e. in the high-temperature region, 
we cannot assume that the coupling between the symmetry and 
the topological defects would be identical as for T ~ T c . Indeed 
with a lot of small Zi walls dividing the system it is not certain 
that the vortex and the antivortex should play an important role 
and we could get a standard Ising transition. 

This hypoth esis w as advanced by Garel et Doniach in 1980 for the 
J1-J2 modelPS 

but numerical MC simulations have shown that 
it is not the cas for this system. 

To choose the most probable hypothesis we will look now at the numer- 
ical results: 

Some information can be provided by numerical MC simulations in the 
high-temperature region. Fitting the results it is possible to get the critical 
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temperature and the critical exponents. Because of the presence of correc- 
tions and by consequence of many free parameters, the extraction of the 
results is difficult. 

For the Ising symmetry Olsson^SI nas found a standard ferromagnetic 
one with v Ismg ~ 1) while Jose and RamireJCSl f oun d ^Ising _ .87. New 
simulations should be done to resolve this discrepancy. If v Istng ^ 1, it will 
a counter-proof to Olsson's claim (first possibility). On the other hand, if 
yising — 2, we cannot rule out the second possibility. 

The expo nent of the KT transition have been only calculated by Jose 
and Ramirez ! ^ ^ They found non standard v KT ~ 0.3 to be compared to 
0.5 in the ferromagnetic case. If this result is in favor of the second hypoth- 
esis, it cannot rule out the first one. Indeed we can reverse the argument 
given at the end of the second hypothesis. If at high temperatures the sys- 
tem is composed of many domains separated by walls of Z%, the correlation 
between the topological defects could not be the one of the standard one. 
Therefore there is no definitive conclusion from the high-temperature sim- 
ulations. 

Now we can look at the results in the FSS. 
For the Ising transition the cr itical ex ponents have been precisely de- 
termined by Loison and Simon , 1 11 for the J1-J2 model using both 
the FSS and the dynamical properties including the first correctio nal^ 
due to the lattice size L. They obtained non standard critical exponents 
yismg _ 0.815(20). The "large" error is due to the inclusion of the cor- 
rections. See Ref. [Ill] for a review of all the numerical works. For the 
triangular lattice the critical temperature is known with high precision 
T Ising = 0.5122(1). 

The KT transition is more problematic. First we can use the helicitjdiniT 
to find the critical temperature. If we admit the same universal jump at the 
critical temperature as the ferromagnetic one, we get T^ T — 0.5010(10) a 
temperature much smaller than T c 7sms = 0.5122(1). But we cannot be sure 

j-i -1- t frustrated ^ f errornaqnetic 1 ,1 • 1, < 1 , 1 .,i 

that I j ump = 1 j ump and this result must be taken with care. 

As explained in the previous section we can use the Binder parameter or 
the dynamical properties of the model to calculate T^ T and the critical 
exponents. Using this method we found very interesting results. First the 
Binder parameter shows a power-law transition and not an exponential 
one. In addition the critical exponents and temperatures calculated by 
the two methods are in good agreement with T^ T = 0.5102(2) just be- 
low Tl smg — 0.5122(1) for the triangular lattice. The critical exponent 
j] KT = 0.36(1) is very different from the standard one 0.25. These results 
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are in accord with those on the J\-Ji model. What is disturbing is that 
for the triangular lattice, accepting T^ T = 0.5102(2), the helicity jump is 
smaller than the Y^[° ma9 " e * IC although it is believed that it should be 
greater. 

The situation is puzzling with no definitive conclusions. 

5.4. Frustrated XY spin systems: Z3 (g) SO(2) 

Two-dimensional systems with an order parameter like Z3 ® 5*0(2) are less 
numerous than three-dimensional systems. Following Table |2 we know at 
least two systems: the triangular lattice with intermediate NNN interaction 
and the Potts- V2,i model. 

No num erical studies have been done for the former system, but we have 
studiecPH the latter one. It gives a first-order transition. 

Following the reasoning of the previous section two principal hypotheses 
appear: 

1. Olsson hypothesis: KT transition at lower temperature following 
the transition associated to the Potts symmetry Z3 at higher tem- 
perature. Non standard critical exponents should appear in the 
FSS region due to a screening length {£ KT ). The "true" stan- 
dard second-order behavior will appear only for very large lattices 
(L ^> S, KT ) or in the high-temperature region. 

2. A new behavior for the KT and Potts transitions at the same crit- 
ical temperature. 

The results are in favor of the second hypothesis. The Potts symmetry 
has a first-order transition. It means that its correlation length £z 3 is finite 
at the critical temperature and less than S;kt- There will be no change 
even if the system size L is much larger and the system will never show 
a standard three-state Potts second-order transition. Obviously we cannot 
apply directly this result to the Ising-KT transition since the two models are 
not equivalent although it is nevertheless an argument against the Olsson's 
hypothesis. 

5.5. Frustrated XY spin systems: Z 2 g) ® SO(2) and 
Z 3 ®Z 2 ® SO (2) 

Following Table [2 the order parameter Z2 <£> Z2 <8> 50(2) can appear with 
a J1-J2-J3 lattice or an Ising-T/2,2 model. In addition Z 3 ® Z 2 ® 50(2) 
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exists in the triangular antiferromagnetic with NNN interaction and in the 
Potts-V2,2 model. No numerical studies have been done on these models. 

5.6. Frustrated Heisenberg spin systems: SO(3) 

For a non frustrated Heisenberg spin system the order parameter will be 
SO(3)/SO(2) and no topological defects exist in two dimensions, no phase 
transition will appear. 

For a frustrated system many order parameters can exist (see Ta- 
ble EJ, but only the planar GS have received attention. In this case the 
order parameter is SO(3) and there exist point defects that is ^-vortex- 
antivorte:;ff2EI2l These topological defects are different from the Z-vortex 
present in XY systems. The existence of a critical transition driven by the 
unboundi ng of vortex-antivortex was first conjectured by Kawamura and 
MiyashitaPH 

At low temperatures the spin waves will dominate the behavior of the 
system and forbid an infinite correlation length below T c contrary to XY 
spins. The behavior should be e quiv alent to the one present in the four- 
dimensional ferromagnetic system!^ This conjecture based on the non lin- 
ear a model (d = 2 + e expansion) and on symmetry arguments was checked 
by Southern and YounJ^ and Caffarel et al.ESI 

At highe r temp eratures the topological defects will clearly have a 
rol J14MMl4SllIfil 

Caffarel et ali ^^ have studied numerically two models having identical spin 
waves but not the same topological defects. One has the topological Z2 de- 
fects, the other not. They showed that the two models are equivalent at 
low temperatures, but show differences at higher temperatures due to the 
topological defec ts. 

Southern and XvJEHI studying by MC simulations the vorticity associated 
to the vortex-antivortex proposed that the vorticity has a jump at the crit- 
ical temperature simil ar to that of the KT transition. 
Kawamura and Kikuch^^l using different boundary conditions in MC sim- 
ulations also observed various phenomena associated to the vorticity which, 
according to the m, ar e a proof of a phase transition. 

Last, Wintel et alP^SIhave studied theoretically and numerically the region 
above the suspected critical temperature for the triangular antiferromag- 
netic lattice (T c ~ 0.29). They claimed that the correlation length and 
susceptibility must follow a KT law. 

There are numerical evidences of the importance of topological defects 
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at finite temperatures and it seems that these systems undergo a kind a 
topological phase transition. Nevertheless we have no certainties that the 
phenomena present at finite sizes would hold for the infinite size. 

5.7. Frustrated Heisenberg spin systems: Z 2 ® SO(3), 
Z 3 ® SO(3) . . . 

By inspecting Table 03 it is not difficult to see that other order parameters 
exist in frustrated systems. One of the most interesting should be Z 2 <£> 
SO(3). It corresponds to a non planar GS which could exist in experimental 
systems. Numerically the simplest system would be the Stiefel V3.3 model 
or equivalently the Ising-V^ model. The comparison with the XY case 
would be very interesting. In particular the Ising transition could appear 
near the SO(3) transition (see previous section) and the coupling between 
the Z 2 walls and the vortex could be very instructive. If we follow the second 
hypothesis of Olsson for the XY case (see section l5T31l , the transition should 
display an "almost second-order transition" . 

Symmetrically the coupling of the SO(3) vortex with a Z 3 Potts model 
which appear in the Potts-V^ model could also disclose interesting prop- 
erties. Does it show first-order properties as for 5*0(2) vortex? 

Some other breakdowns of symmetry could also appear like Z$ ® Z2 ® 
SO (3) in the antiferromagnetic triangular lattice with large NNN interac- 
tion. Similar questions appear. We note that the SO(3)/SO(2) order pa- 
rameter does not have a topological transition. The transition of the type 
Z 2 <g> SO{3)/SO{2) for the J x -J 2 model or Z 3 ® SO(3)/SO(2) for the tri- 
angular lattice with intermediate NNN interaction should have a standard 
transition (Ising or three-state Potts model). 

5.8. Topological defects for N > 4 

Some other topological defects should exist for an order parameter of the 
type SO{N), with a GS in N - 1 dimensions, Z 2 <g> SO{N), with a GS in TV 
dimensions, and Z q (g> SO(N) with a coupling of a GS in N — 1 dimensions 
and a lattice symmetry • • • All questions raised for Heisenberg spins still 
hold in these cases. 

6. General conclusions 

In this chapter we have studied the phase transition in frustrated systems 
between two and four dimensions. We have found various breakdowns of 
symmetry, contrary to the unique one for ferromagnetic systems. 
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In three dimensions the transition is always of first order in the ther- 
modynamic limit. However for "small" sizes in numerical simulations or 
for temperatures not "too close" to the transition temperatures in ex- 
periments, the system could display an "almost universality class" for an 
0(N)/0(N — 2) breakdown of symmetry. Many compounds studied exper- 
imentally are in this class. 

In two dimensions the situation is much less clear. Indeed the topological 
defects can play a fundamental role and their couplings with a discrete 
symmetry (Ising or Potts) is unknown. We hope that in the near future 
the two-dimensional case will be clarified as the three-dimensional one on 
which our understanding has increased considerably in the last decade. 
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Appendix A: Monte Carlo Simulation 

The purpose of this section is just to show the fastest algorithm to run 
Monte Carlo simulations, the best method to distinguish a first-order tran- 
sition from a second-order one and to extract the critical exponents with a 
reliable estimate of errors. 

Markov chains and algorithms: Since we cannot enumerate all the 
spin configurations, we are forced to use numerical simulations to calculate 
the physical quantities. The method of choice is the Monte Carlo method 
generating a set of phase-space configurations from a Markov chain de- 
fined by the transition probability W[s,s'] between states {s} and {s'}, 
where the new configuration depends only on the pre cedin g one. There are 
many ways to get the transition probability W[s, s']P-^ The Metropolis 
algorithm is the simplest but not ver y effi cient one. A new heat-bath algo- 
rithm, called Fast Linear Algorithm) " I for vector spins is more efficient. 
For example it is three times faster than the Metropolis algorithm for a 
two-dimensional triangular antiferromagnetic latti ce at the critical temper- 
ature. Furthermore the use o f the over-relaxation^^ can reduce strongly 
the autocorrelation time.H^ The cluster algorithm cannot be used for frus- 
trated spin systems. Indeed there are two problems. The first is to take into 
account the competition of different interactions in a plaquette. It means 
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that many spins should be considered in one step when constructing the 
cluster. Second, even if we were able to construct a cluster, the BS is no 
more 0(N)/0(N - 1), but 0(N)/0(N - 2), for example. Cons equently a 
symmetry to a N — 1 plane used in the Wolff's algorithm^ 

should be 

changed. 

Thermodynamic averages are taken simply by 

Nmc 

A = (1/Nmc) £ A[t] (9) 
t—i 

where Nmc is the number of new configurations generated (Monte Carlo 
Steps) and A[t] the value of the quantity at step t. The great advantage of 
this method is that the partition function Z = e.~@ E w needs not be 

calculated. 

To estimate the errors of the averaged quantities we have to take into 
account that each new configuration is correlated with the previous one. 
We define the autocorrelation time r by the number of MC steps required 
to obtain two uncorrelated spin configurations. It is calculated using the 
autocorrelation functiorffS c(t) = (l/ X )[(A(t)A(0)) - (A) 2 ] where time 
is measured in MC steps (see comments in the appendix of Ref. [22]). If 
-/V 3> t, then a useful approximation for the error in A is given by 

^ = AW(l + 2r) < 10 > 
X = (A 2 )-(Af. (11) 

This expression is identical to the standard deviation but with an effective 
number of independent measurements given by 

Nm^ (12) 
1 + 2t V ' 

Problems arise when quantities are a combination of different averages, 

for example the susceptibility (jTT|) . We could try to treat (A 2 ) and (A) 2 as 

independent quantities and estimate the error by the sum of the errors of the 

two quantities. But the result will be overestimated due to the correlation 

between the two elements of the s um. T o solve this problem we can use, for 

example, the Jackknife procedure! ^^ An application of this method can 

be found in the appendix of Ref. [22] . 

The histogram method: The great advantage of the histogram 
method in the analysis of MC data is that a run at a single tempera- 
ture T\ can be used to extract results for a continuous range of nearby 
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temperatures JED i n practice the range of temperature over which (A) may 
be estimated from a single MC run at T\ is limited by the range over which 
reliable statistics can be expected for H(E), the histogram of the energy. 
A rough guide is T a < T < T),, where T a and Ti, correspond to the energies 
(E a ) and (E b ) at which H(E) ~ \H max . Since H(E) becomes more sharply 
peaked with system size, the valid temperature range becomes smaller as L 
increases. Multiple histograms made at a number of nearby temperatures 
may be combined to increase accuracy. 

Nature of the transition: Differentiating a weak first-order from a 
second-order transition could be difficult. The finit e size scal ing (FSS) for a 
first-order transition has been extensively studied j te>^ | too | lo4 |^ first-order 
transition should be identified by the following properties: 

a) The histogram of the energy, Pt(E), has a double peak. 

b) The minimum of the fourth order energy cumulant W varies as: 
W — W* + bL~ d where W* is different from 2/3. 

c) The temperatures T(L) at which the specific heat C or the suscep- 
tibility x has a maximum should vary as T(L) = T c + aL~ d . 

d) The maximum of C and \ are proportional to the volume L d . 

Obviously a double peak of Pt{E) which becomes more pronounced when 
the size increases is preferable (way a). If the two peaks are too close, 
the three other possibilities (b, c and d) can check if the probability is 
gaussian (second-order transition) or not (first-order transition) in the limit 
of infinite size. 

Second-order phase transition: For a second-order phase transition 
the interesting parameters are the critical temperature and the critical ex- 
ponents. There are at least three main ways to calculate them: 

a) Consider the high and low temperature regions where the correla- 
tion length There, we can fit £ = a^(T — T c ) _l/ + corrections. 
The corrections become less important near T c , but then £ is very 
large, and a very big system size is necessary. Furthermore the au- 
tocorrelation time r ~ (T — T c )~ z becomes also very large and the 
e rror 5 £ becomes very big even if we can use some tricks to diminish 

b) The Finite Size Scaling (FSS) region is the most powerful method. 
It is the region with £ (theoretical) 3> L. The best method is to cal- 
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culate the Binder cumulanl -*-^^ ! Ul (T) = 1 — 3 |^ 2 | 2 for different 
size L. 

1. We can calculate the critical exponents directly by plotting 
yjj-ilv as a function of Ul(T) and all the curves must collapse 
for the correct exponent —7/^ = ij — 2 .H53 Indeed we can write 
xL ??_2 = h{Ui,) with h an unknown function. The other exponents 
P/p and v can be obtained similarly using the magnetization and 
^1 = ^\m) ~~ (E)- It is the fastest way to get the critical ex- 
ponents. This method works well even for the Kosterlitz-Thouless 
transition! ^ 

2. We can also use Ul(T) to calculate the critical temperature T c 
using the crossing of Ul for different sizes. Very good statistics is 
needed because the evaluation of the first correction is necessary 
to get a correct value of T c . 

3. Having determined T c we can calculate the critical exponents 
using, for example, \ cx L+t/" , • • • 

4. The last properties are also valid for the maximum (or minimum) 
of x and V\ , but we need many simulations at various temperatures 
to find the maxima because their locations vary as a function of 
the size L and are differe nt for each q uantity. 

c) The dynamical propertie J 1 59 1 1 W | l 6 1| 1 1>2 | 1 fofl are no t very often 
used even if it is surely the fastest method available. One has to 
prepare a state in the GS (i.e. T — 0) or randomly (i.e. T — 00), 
and to observe the dynamical properties of the system at the crit- 
ical temperature T = T c for a finite number of Monte Carlo steps 
before the equilibrium is reached. For an example see Ref. [111]. 

Kosterlitz-Thouless (KT) transition: A KT transition 11^ 1 exists 
for two-dimensional lattices with XY spins. The unbounding of vortex- 
antivortex pairs appear at the critical temperature, T c ~ 0.9 for square 
lattice. 

This transition has some special features: for T < T c the correlation length 
(£) is infinite, while for T > T c it has an exponential decreasing contrary 
to the power law behavior in the standard transition. 
To find T c and the critical exponents we have several ways: 

a) we can fit £ in the high temperature region, which is problematic 
because of the many free parameters and of the exponential form. 

b) we can use a method using the behavior of various quantities in 
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the finite size scaling region (FSS) where £ ^> L. 

1. The first one is to use the universal jump of the helicity T at the 
critical temperature (in the finite size scaling region where £ 3> £), 
but this method requires the jump of T(L) which is known for the 
non frustrated case, but not for the frustrated case. 

2. It is therefore interesting to find another method to get the 
critical exponent, even without calculating T c . The method (b.l.) 
introduced previously for a second-order phase tra nsitio n using the 
Binder parameter works well for the KT transitionPSSlThe Binder 
parameter, contrary to the common belief, crossed around T c . Plot- 
ting xL n ~ 2 as a function of Ul for various sizes, the curves collapse 
at the correct value of r\ = 0.25. This method can be applied what- 
ever the form of £ is. Therefore it should be applicable to all phase 
transitions. 

c) The dynamical behavior of this system can also be used 

For the two-dimensional ferromagnetic XY spin system, all methods men- 
tioned above are in agreement. 

Appendix B: Renormalization Group: Landau-Ginzburg 
theory, expansions in fixed dimension d — 3 and for 
d — 4 — e and its implications for experiments 

In this short section we would like to give the necessary knowledge to help 
to understand the concepts of fixed points, flow diagrams and crossovers. 
First consider a system of Heisenberg spins with a Hamiltonian: 



and the spin is restricted to be of norm 5 = 1. Replacing this constraint by 
an exponential potential S = 1 f e u ^ s ■ dS the Hamiltonian can 
be written in the form 



H = KY / (S(i) - S(j)) 2 + r J2 S(*) 2 +"o£ (S« 2 ) 2 . (14) 



Two remarks: 

1. Other terms could be added in the Hamiltonian. These additional terms 
could be unimportant that means they become irrelevant near the fixed 
point (see later). Or we cannot treat them, that means they are not renor- 
malizable in technical terms. If these neglected terms are important the 




(13) 



<ij> 
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method cannot describe the physics of the phase transition. This unfor- 
tunate case occurs for frustrated system for low dimensions of space and 
spins 

2. In addition, since the transformation is not exact, the starting value 
cannot be known even if it is possible to make a rough guess 
We have plotted in the Fig.^lthe RG flow for this model. 







, T 

T c (u) 
















slow 





Fig. 13. RG flow for 1141 . F is the stable fixed point, G is the Gaussian one. Near F 
the flow is very slow. 

The system will physically keep its value uq. But its RG flow will not. 
Following a series of steps the value of u will change in approaching the fixed 
point F. At each step the correlation length of the system will be divided 
by a factor b > 1. In addition the flow becomes slower near the fixed point 
F. Consequently, to reach the neighborhood of F, where the system have a 
second-order transition, the initial £ must be large "enough" . We know that 
£(T) ~ (T — T c )~ v for a second-order phase transition. Hence, to obtain 
numerically or experimentally £ which follows a power law (without too 
many corrections), the temperature must be close "enough" to the critical 
temperature T c . Numerically, in the Finite Size Scaling (FSS) regions where 
£ (theoretical is much bigger than the size of the system L, the FSS law (see 
EJ) will be valid only for L large "enough" . We observe that in the Fig. 
only the plane T c (u,v) is plotted. 

Consider now a model with a XY anisotropy, for example we add to 
the Hamiltonian (O a term DY Jl {St) 2 with < D < 1. We have now 
two parameters in our Landau-Ginzburg model: u for the length of the spin 
and v associated to D. The critical plane (u, v) of the flow diagram of this 
model is plotted in Fig.^J From the initial point (uq, vq 1) the flow will 
go close to the Heisenberg fixed point Fjj and then has a crossover to the 
Fxy fixed point. Near Fh the flow is very slow and needs a lot of steps 
to escape from the influence of Fh and reaches finally the neighborhood of 
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Fxy- Therefore, to observe the "true" XY behavior, the correlation length 
must be very large. That means that the temperature must be very close 
of T c or, if we are in the FSS region, the system size must be very large. 




Fig. 14. Schematic RG flow (T c (u, v) plane) for ITHwith O^^S, 2 ) 2 ). F H is the Hciscn- 
bcrg fixed point, Fxy is the XY fixed point, G is the Gaussian one. Near Fjj and Fxy 
the flow is very slow. 



We must know the stability of the fixed point when anisotropics are 
present. The "quasi" general rule is that the system goes to the fixed 
point with the lowest spin symmetry and the biggest space dimensions. For 
example consider the quasi one-dimensional system CsMnBr3. This com- 
pound will display a one-dimensional behavior before reaching the three- 
dimensional behavior. Another example is Cu(HCOO)22CO(ND2)22D20 
which has a small Ising anisotropy. In this case the system will show a 
Heisenberg behavior before showing the Ising behavior. In real compounds 
small "anisotropics" are almost always present. Therefore many crossovers 
will appear when the temperature approaches the critical temperature. 
Consequently interpretations of experiments could be difficult. 

To get the critical exponent and the flow diagram, the expansion in 
fixed dimension d or in d = 4 — e consists of expanding the exponential 
e Hamiitoman{u) arounc j Gaussian fixed point u — 0. Results are in the 
form of a power series of u. Unfortunately this series is not convergent, but 
in the ferromagnetic case it has been shown to be resummable with Pade- 
Borel techniques. However, there are many ways to resum and one must 
be chosen following some criteria!^ This resummation has been proved 
efficient for ferromagnetic systems, however it is not certain that it will work 
for frustrated systems. Indeed the presence of two vectors (O(N) /0(N — 2) 
BS in this case (compared to the 0(N)/0(N — 1) BS in the ferromagnetic 
case) gives rise to a series of power of u and v (u,v, uv, u 2 , • • • ) and this 
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double expansion is difficult to resum. 
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